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I.   INTRODUCTION 

The  various  electronic  structure  models  used  by- 
chemists  to  help  them  understand  and  describe  chemical 
behavior  have  changed  very  little  in  the  more  than  fifty 
years  since  the  formulation  of  the  first  principles  of 
quantum  mechanics.   The  concept  of  valence  electrons  being 
responsible  for  most  of  the  chemical  behavior  that  has 
been  studied,  is  one  which  has  always  stayed  alive  and  is 
now  again  coming  into  the  foreground.   Hov/  the  Is  orbital 
of  a  transition  metal  changes  when  forming,  say,  a  hydride 
is  not  a  question  which  does,  nor  need  concern  us  since  we 
all  know  this  change  is  slight.   Our  problem  is  that  quan- 
tum mechanics  and  the  principle  of  indistinguishable  elec- 
trons has  placed  upon  us  the  burden  of  having  to  include 
core  electrons  in  our  calculations  in  order  to  properly 
account  for  the  electronic  structure  of  atoms  and  molecules. 
In  order  to  carry  out  calculations  on  molecular  systems, 
using  Ilartree-Fock-Roothaan  type  equations,  various  highly 
sophisticated  computer  programs  have  been  developed  and 
refined  to  stages  where  they  can  handle  virtually  any  mole- 
cular system  of  interest.   As  more  complicated  systems  are 
studied,  the  computer  tir,v_  ..;.u  .;pace  required  to  perform 
these  calculations  places  an  upper  limit  on  our  capabilities 
To  overcome  this  problem  various  approximation  techniques 
have  arisen.   Those  techniques  all  have  a  common  theme:   it 


is  the  valence  electrons  which  undergo  significant  changes 
during  reactions.   Unfortunately  it  is  generally  found 
that  the  simpler  the  calculations  become,  the  greater  the 
loss  in  accuracy  becomes.   The  present  v;ork  contains  a 
method  which  also  leans  on  the  concept  of  valence  electrons, 
and  results  in  a  great  reduction  in  the  size  and  cost  of 
calculations,  but  which  is  found  not  to  be  accompanied  by  a 
3erious  loss  in  accuracy. 


II.   BACKGROUND 

In  order  to  appreciate  the  method  used  here,  a  brief 
account  of  earlier  works  v/hich  led  to  it  is  presented. 
This  review  is  not  intended  to  be  an  historical  one,  nor 
a  complete  survey  of  all  the  v/ork  done  to  the  present  on 
pseudopotentials,  effective  potentials,  model  potentials, 
etc.   The  intention  here  is  to  present  some  of  the  results 
of  earlier  investigations  into  approximation  techniques 
which  in  some  sense  led  to  the  present  scheme  as  both  a 
logical  and  intuitive  conclusion. 

A.   Pseudopotentials 

From  experimental  and  theoretical  studies  on  numerous 
molecular  systems,  we  know  that  it  is  the  valence  orbitals 
v;hich  are  responsible  for  chemical  bonding  etc.  ,  and  that 
it  is  primarily  the  outermost  part  of  the  valence  orbital 
in  which  this  takes  place.   It  is  therefore  important  that 
any  method  developed  to  predict  chemical  properties  give 
an  accurate  description  of  the  "tail"  of  the  valence 
orbitals.   The  other  region  of  the  valence  orbital,  the 
"core"  part,  also  requires  careful  consideration,  since 
this  is  the  part  that  must  be  kept  out  of  (orthogonal  to) 
the  core  orbitals  in  order  that  we  arrive  at  physically 
realistic  results.   These  two  features  of  valence  orbitals, 
although  of  the  utmost  importance,  have  an  important  dif- 


ference.   For  the  tail  region  of  the  orbital  it  is  not  only 
important  that  this  part  be  accurately  accounted  for,  but 
that  we  be  able  to  retrieve  the  information  that  describes 
that  region.   On  the  other  hand,  for  the  "core"  region  of 
the  valence  orbital,  we  only  require  that  it  be  handled 
with  sufficient  accuracy  to  prevent  physically  unrealistic 
results.   This  has  been  taken  advantage  of  by  Weeks,  Hazi 

and  Rice   who  have  given  a  generalization  of  a  method  due 

2 

to  Phillips  and  Kleinman   known  as  the  pseudopotentxal 

method.   Since  the  core  region  of  the  valence  orbital  re- 
quires extra  basis  functions  for  the  reproduction  of  nodes 
and  wiggles;  and  since  there  is  no  real  interest  in  having 
an  accurate  description  of  these  nodes  and  wiggles,  their 
only  purpose  being  to  maintain  orthogonality  to  the  core, 
their  method  allows  one  to  use  a  smooth  pseudowavefunction 
for  the  valence  orbital.   Since  it  is  a  smooth  wave function 
it  requires  fewer  basis  functions  for  its  description,  lead- 
ing therefore  to  a  simplification  of  the  calculations. 
Mathematically  they  transform  the  Hamiltonian 


Hd)   =  (T+V)({)   =  e  (})  (la) 

^v        ^v    v^v 


subject  to 


<4)  U  >  =  5  (lb) 

^c'^v     cv 


to        H^^x   =  (T+V+V^^'^'^^')X   =  e  X  (2) 

Ay  / Ay  V   V 


v/here 


^PS,GPK  ^  _^jp  _  pj,  _^  p^p  +  e  P  (3) 


^  -    ---    --    -V 


is  the  "Generalized  Phillips-Kleinman"  pseudopotential, 

and 


core      core 
P  =   I   P^,  =   I    Icxcl  (4) 

c        c 


is  the  projection  operator  v/hich  projects  onto  the  core 

PS 
orbitals.   The  Phillips-Kleinman  form  for  V  -  is 


PS,PK    ^°r^®  ,         ,^  ,c\ 

V      =   )  (e   -  e  )P   .  (5) 

^     V    c   c 
c 


The  pseudoorbital  above,  x^,    need  not  be  identical 
to  4)  .   c   will  be  an  upper  bound  to  e  for  any  arbitrary- 
basis  set.   Thus,  we  are  at  liberty  to  choose  for  x^  ^o^® 
small  basis  set  which  contains  only  those  functions  re- 
quired for  an  accurate  description  of  the  tail  region.   To 
see  how  e   is  an  upper  bound  to  c,  we  note  that  any  arbi- 
trary function  x   can  be  v/ritten  in  terms  of  the  complete 
set 

Y=ya*+<})+ya't'  (G) 

''v    '•   c  c    V    ^   e  e 


Any  component  of  Xw  ^^^   *^°  ^^^   first  term,  has  its  energy 

lifted  to  c   by  the  pseudopotential  which  adds  an  .; 

c   -  e   onto  E   (for  6  )  giving  e..  as  it..  v...^igy.   The 
V     c        c        C  V 


second  term  already  has  e   as  its  energy  and  any  conponent 
of  y   from  the  third  term  will  be  higher  than  e   so  that 
any  arbitrary  x  will  in  fact  have  e   as  an  upper  bound. 

It  should  be  noted  that  in  the  Hainiltonians  described 
above,  no  approximations  have  been  made,  but  rather  a 
transformation  from  an  equation  (la)  subject  to  constraint 
(lb)  to  a  new  unconstrained  equation  (2) .   These  are  the 
only  Hainiltonians  correctly  termed  "Pseudohamiltonians" 
whereas  those  approximations  to  be  presented  below  are 
"model"  or  "effective"  Hamiltonians. 

The  problems  with  the  above  Hamiltonians  are  clear. 
First,  if  care  is  not  taken  it  is  possible  to  choose  a 
X   which  has  very  little  resemblance  to  <j)   except  in  energy, 
Second,  consider  the  ground  state  of  Na.   In  order  to  solve 
for  the  3s  orbital  above,  we  must  know  the  core  orbitals 
which  we  must  get  by  a  Self  Consistent  Field  type  calcula- 
tion which  would  have  given  us  the  3s  orbital  any\>;ay.   To 
solve  for  an  excited  state  we  must  either  recalculate  the 
projection  operators  P   or  assume  they  have  not  changed, 
whereby  the  method  becomes  tantamount  to  a  frozen  core  ap- 
proxim.ation.   In  addition,  v;e  are  left  with  the  undesire- 
able  coulomb  and  exchange  integrals  over  the  core  orbitals. 

B.   Semiempirical  Effective  and  Model  Potentials 

To  help  improve  the  situation  a  number  of  approxima- 
tion techniques  have  been  developed.   For  convenience  Table 
I  gives  some  of  their  Hamiltonians  along  with  the  Phillips- 


Kleinman  Pseudohamiltonian.   Note  that  although  some  of 
these  have  been  developed  for  many  valence  electron  systems, 
only  the  one  valence  electron  Hamiltonians  are  presented, 
allowing  for  ease  of  comparison.   For  more  than  one  val- 
ence electron,  the  appropriate  coulomb  and  exchange  opera- 
tors over  the  valence  orbitals  must  be  added. 

The  first  class  of  these,  the  semierapirical  poten- 
tials, are  usually  formed  by  adjusting  the  parameters  to 
fit  experimental  binding  energies.   Although  this  corrects 
somevvhat  for  correlation  effects,  there  is  no  guarantee 
that  accurate  wavefunctions  and  other  properties  will  be 
correctly  reproduced. 

Sim.ons  and  Mazziotti   have  tried  using  a  Gaussian 
basis  set  v/ith  the  Hamiltonians 


H  =  T  -  ^  -f  ^  e-2^^  (7) 

r    r 


and 


2 

H  =  T  +  C,r   +  C^  (8) 


and  found  them  to  give  poor  results  for  molecules  due  to  a 
lack  of  symmetry  (S.)  dependence.  In  addition,  C-,  of  equa- 
tion (8)  does  not  die  off  as  it  should  on  say  another  atom. 

7 
The  Abarenkov  and  Heine   potential  improves  on  this  some- 
what by  introducing  an  i.   dependence  and  a  cutoff  point,  i.e, 


H  =  T  +   );   B„P,    r  <  R  (9a) 


TABLE  I.   One  Valence  Electron  Hamiltonians 


Pseudohamiltonians 


Examples 


n   core 

V 


H  =  T ^+  y   (G,-e^) 

r    '-     V   c 


Phillips-Kleinman' 


Semi-Empirical  Effective  Kamiltonian 


H   =   T   -   ^  f  ^  e-2^^ 
r         r 


3  4 

Hellmann,      Gombas, 

Szasz    and  McGinn, 


Simons    and  Mazziotti 


.6 


H   =    T    +   C^r      +    C2 


Simons   and  Mazziotti 


H   =   T   +      I      Pjj^B£    r<R      =   core   radius 


T  - 


n 

V 

r 
n 


H  =  T 


H  =  T 


r    ^ 


r>R 


^^£ 


n 


£   r 
A, 


_  ^  +  y  -^ 
r    ^   r 


e   ^    P, 


Abarenkov  and  Heine 


Simons  ,  Goodfriend 

9 
and  Hart 


Schwartz  and  Switalski 


.11 


n 

H  =  T  -  ^ 


A  -ar     atoms , 

—  e 

r      a  >  S,max  of  core 


£max  of  core 
£  =  0 


P„B„  e 


-6  r 


atoms , 
I    <   ilmax  of  core 


Jimax  of  core 
Z    =    0 


PnC„  e 


■^£^ 


molecules 


Chang  et  al 


12 


TABLE  I   Continued 


Ab  Initio  Ef f ective-Haniltonians 


Examples 


H  =  T  -  "^  -  ^  e""^  +  I  (e  -e  )P 
"        r    r         '-vcc 

c 


n 

=  T  -  —  +y  (e-e:)P   No  6  on  sane 

r    ^  ^  V  c   c     I 

c  center  as  v 


Coffey  et  al 
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n   n  2 

n    a   V     -  a  _.  r 


H  =  T 


r  .   r   -'-  '•   c  c 


H  =  T  - 


X 

c 

n 

V 

r 

+ 

ilval 

I 
1  =  0 

I 
i 

^t'^i, 

i/' 

n . 

1 

Bonifacio  and 
Huzmaga 


Durand  and  Barthelat 


15 


H  =  T  - 


n     £val 
^    £  =  0  i     '^ 


'   e 


Kahn  and  Goddard 
Melius 


16 


n 

H   =   T-—        r>R  (9b) 

r  c 


where   P-    is   the   angular  part  of   the   projection  operator 

P^   =    I    Umxilml    =    I    I^j,j„><Y^^U  (10) 

m  m 


and  where  R   is  the  radius  of  the  core  orbitals,  and  n 

is  the  number  of  valence  electrons.   The  Simons   potential 


—  +      I     -^  (11) 

^    £=0  r 


dies  off  smoothly  without  requiring  a  cutoff,  and  contains 
an  I  dependence  of  the  valence  orbitals.  There  is  also  an 
analytic  expression  for  the  energy 


-2z2 

-  =  0,1,2,. ..  .    (12) 


^         {2P+1[(22+1)^  +  8B^]-'-/^}^ 

9 

Goodfriend  and  Hart   have  developed  some  of  these  poten- 
tials by  setting  P  =  0  and  using  as  a  wavef unction: 

1/2 


>:v 


e^ 


3tt 


re  ^^  (13) 


where  £  \\?as  obtained  by  fitting  to  the  STO '  s  of  Clementi 
et  al_.  To  get  the  nth  excited  state,  P  was  then  set 

equal  to  n. 

Schwartz  and  Switalski    used  what  they  call  a 
"Gaussian  Screened  Coulomb  Potential" 


10 


2 
n         A.   -Y„r 

^  +  I  P,  ^  e  '  (14) 


which  introduces  an  J,  dependence  in  a  Gaussian  term  which 
allows  for  easier  use  in  molecular  calculations  and  in- 
corporation into  existing  programs.   They  have  also  tested 
the  effect  of  taking  the  original  basis  set  from  which 
the  potential  was  developed  and  leaving  out  the  first  few 

core  type  functions. 

12 
Chang  et  a_l.    have  also  investigated  a  screened 

coulomb  potential  with  a  basic  form  given  by 


V  _  A  p-ar 
r    r 


For  the  first  row  they  find  the  potential  to  be  fairly 
insensitive  to  the  value  of  the  parameter  A.   As  they  in- 
crease the  atomic  number  the  results  get  poorer.   For  Ti 
for  example  they  find  that  they  require  a  polynomial  ex- 
pression for  A.   They  attribute  this  to  the  greater  error 
in  the  frozen  core  approximation  for  larger  atoms. 

C.   Ab  Initio  Effective  and  Model  Potentials 

The  second  class  of  effective  Hamiltonians  consist 
of  the  theoretically  based  Hamiltonians,  that  is,  those 
which  do  not  make  use  of  experimental  values  for  parameter 
adjustment.   The  simplest  one  in  Table  I,  that  of  Coffey 
et  al^.  , 

11 


H  =  ■ 


n    n   _ 

T ^  e  "^  +  y  (e  -e  )P               (15a) 

r    r         ^  V  c   c 
c 

for  both  orbitals  on 

^v    Y  /      \  TD  ^^  atom  different     m  c:k\ 

T  -  —  +  Z  (e^  e^)P^  ^^^^   ^^^   potentials  ^^^^^ 


consists  essentially  of  the  Phillips-Kleininan  potential 
•yith  a  parameterized  exponential  term  to  replace  the  cou- 
lomb and  exchange  terms  of  the  core  orbitals.   They  have 
found  this  form  to  be  basis  set  dependent,  i.e.,  it  must 
be  readjusted  for  different  basis  sets.   They  also  find 
that  when  core  functions  are  included  one  pseudoorbital 
can  collapse  into  the  core  when  trying  to  avoid  another. 
Thus  the  result  of  a  calculation  with  additional  functions 
in  the  core  region,  which  should  give  better  results 
(greater  flexibility)  gives  poorer  results.   They  also 
find  that  e   (in  the  projection  operator)  must  be  readjusted 

after  each  iteration  in  an  SCF  type  calculation. 

14 
The  Bonifacic-Huzinaga    Hamiltonian 


n      a  n  A.  -a.r     core 

H  =  T  -  --^  -  y   -^^-i  e   ^    +   y   B  P         (16) 

r    V    r  i-    c  c 

1  c 


uses  a  series  of  Gaussians  with  adjustable  exponents  and 
coefficients,  and  changes  the  e  -e   part  of  the  Phillips- 
Kleinman  type  term  to  B   (an  adjustable  parameter) .   For 

first  row  atoms,  however,  the  B   parameter  matches  c  -g 

14 
very  well    so  that  it  might  have  been  left  at  the  pre- 

12 


dieted  c    -e  .   The  potential  does  have  the  advantage  of  hav- 
ing a  Gaussian  form  which  leads  to  analytic  multicenter 
integral  expressions,  and  can  thus  be  incorporated  into 
existing  Hartree-Fock  programs  (in  fact  the  author  has 
adjusted  the  programs  for  effective  potentials  described 
below  to  run  with  this  type  of  a  potential,  but  has  as 
yet  only  performed  some  preliminary  tests  on  them  for  atoms 
and  small  molecules).   Bonifacio  and  Huzinaga  have  found 
this  potential  to  be  basis  set  dependent.   The  B   parame- 
ter must  be  readjusted  for  different  basis  sets.   They 
also  find  that  for  the  first  rov;  the  11-S  type  Gaussian 
primitive  functions  can  only  be  reduced  to  six,  by  using 
their  potential,  thus  leaving  the  most  important  core  type 
(large  exponent)  Gaussians. 

All  of  the  above  cited  theoretically  based  potentials 
have  their  i   dependence  in  the  core  projection  tern  P  . 
Note,  however,  that  the  effective  potentials  must  in  part 
replace  the  exchange  integrals  between  the  core  and  val- 
ence electrons.   Thus  the  effective  potential  must  contain 
the  ?,  dependence  of  the  valence  orbital.   For  example,  for 
the  first  row  neither  a  2p  nor  a  3d  orbital  has  any  pro- 
jection on  the  Is  core  orbital,  so  that  in  the  above  schemes 
they  would  feel  the  same  potential.   However,  since  the  ex- 
change operator  depends  on  its  operand,  the  2p  exchange  in- 
tegral has  an  i  =  1  dependence  whereas  the  3d  exchange  in- 
tegral has  an  Ji  =  2  dependan  •  . 

15 
Durand  and  Darthelat    have  recently  presented  an 


1? 


effective  Hamiltonian  which  has  the  proper  I   dependence; 


where 


n 

H  =  T ^  +  W_  (16a) 

r     R 


I   n      r 


or  for  many  electrons 


n 

H  =  T ^  +  V,   +  W-,  (17) 

r     4)^    R 


where  V.   are  the  appropriate  coulomb  and  exchange  terms 

over  the  valence  pseudoorbitals.   They  generate  W^  by 
first  finding  the  intersection  of  0   with  the  highest 
core  orbital  of  the  same  symmetry  (R   =  point  of  inter- 
section) then  fitting  for  a  x^  by  minimizing 

<y   -  (1)  !  Y   -  *  >^  (18) 

'^V      V  '   V      V  R   ->■  00 

c 


subject  to 


<X   X  >  =  !• 

V   V 


By  choosing  x   as  a  linear  combination  of  Slater  orbitals 


Xv  =  },    ^i^i 
1=1 
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they  arrive  at  the  system  of  linear  equations 


1 1 

y   [<s.|s.>_-x<s.ls.>]c.  =<S.|y>, 


c 


(see  Appendix  I  for  a  similar  derivation) .   To  obtain 

their  effective  potential  they  must  then  determine  the 

parameters  A     (Eqn.  15b) .   This  is  accomplished  by 
n  f  X. 

.•ninimizing  the  norm  of  the  operator 


O  =  I    [(e:;|x;><x;I)  -  (sfXv><X^I)] 


where  y   is  the  pseudoorbital  from  Eqn.  (18)  above,  e 
is  the  valence  orbital  energy  from  a  complete  SCF  calcu- 
lation and  where  e'  and  x'  are  the  orbital  energy  and 
pseudoorbital  resulting  from  a  particular  set  of  A 

As  a  result  of  their  investigations  of  these  effec- 
tive potentials  for  the  first  and  second  row  of  the  perio- 
dic table,  they  predict  that  to  get  improved  accuracy  they 
will  have  to  replace  equation  (16b)  with  a  more  elaborate 
form. 

The  last  Hamiltonian  in  Table  I  has  the  same  form  as 
will  be  used  in  this  work  (this  v/ill  be  described  separate- 
ly below) .   It  was  first  used  by  Kahn  and  Goddard,    and 
has  the  desired  angular  momentum  dependence  described 
earlier.   For   one  valence  electron  we  have: 


n  rrf 

H  =  T +   U  (19) 

r 


where 

.eff 


U     =  HUo  (r)  UmxJlml  (20a) 

LMAX-1 
=  "lMAX^^^  \l    ,    U^,^^)-\mAX^^))\^^^><^^\     (20b) 


They  obtain  the  potential  by  solving  for  a  Gl  wavefunction 

eff 
(i.e.  X  )  ^ri<^  inverting  equation  (19)  to  find  U    ,  that 

is  they  operate  on  x   with  H  and  then  divide  by  x^  ^^•'^ 

solve  for  U 

,^eff       ,  z  .     V  ,~,  . 

"    =  ^v  +  ?  -^  ^X"  •  ^^^^ 

V 


The  Gl  wavefunction  is  used  because  it  is  nodeless  and 
thus  does  not  lead  to  singularities  during  the  inversion. 
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III.   THE  CORELESS  HARTREE  FOCK  AB  INITIO  EFFECTIVE  POTENTIAL 

The  present  method,  the  Coreless  Hartree  Fock  Effec- 
tive Potential  (CKFEP) ,  developed  by  '  ^rl  F.  Melius  et_ 

17-19 i29 
al.         is  designed  to  include  the  important  features 

rrentioned  above.   This  method  is  now  presented  in  detail. 

A.   The  Coreless  Kartree  Fock  Orbital 

We  v/ish  to  develop  a  potential  which  will  allov/  us 
to  use  smooth  valence  orbitals  (therefore,  fev;er  basis 
functions)  for  our  SCF  calculations.   A  starting  point 
then,  v;ould  be  to  take  such  a  smooth  orbital  and  find  v;hat 
potential  is  needed  to  generate  it.   This  smooth  orbital 
should  be  essentially  identical  to  the  original  valence 
orbital  in  the  tail  region  v/hile  at  the  same  time  be 
"smoothed  out"  in  the  core  region.   This  orbital  v;ill  not 
have  the  core  character  of  the  original  wavefunchion  and 

for  the  orbitals  generated  here,  is  termed  the  "Coreless 

17-19,29 
Hartree  Fock"  (CIIF)  orbital.   (C.  F.  Melius  et  al.)  . 

The  CHF  orbital  is  given  as  -: 


n 

CHF    .HF  ^  v^    aHP  /',^> 

X     -    <i)         +    }       ei    <t>  (22) 

A^      ^v'cc 
c 


which  gives  upon  renormalization. 


n  +1 
Xv    =    I      "i*i  ^23) 
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v;here  c"^  and  6  "   arc  the  core  and  valence  orbitals  re- 
c       ■  V 

spectivelv,  and  n  is  the  n'jraber  of  core  orbitals.  To 
detemine  y  "  '  ci-  ^^  adjusted  to  give  t?ie  best  fit  to 
an  n  =  3  Slater  orbital 


2   -  cr 
K   r   e  ^ 
s 


(24) 


where   IT      is    a  nojcmalization  constant.      This   is   accomplish- 
s 

ed   by  rainirp.izing   the    square   of    the   difference   betv/een.   y' 
and   S ,    i.e. 


3a, 


HF 
/    c;  ■  p  . 

(.1 


2  ^ 

w(r)r      d' 


=    0 


(25) 


V7ith 


tj 


(r) 


m 


(26) 


and  v;here  n  is  originally  set  to  zero  (i.e.  no  v/eight 
factor)  ,  and  C  i^^-  equation  (24)  is  set  at  /-2e    ,    and 
solving  for  >;'  (Hef-erence  29).  We  then  calculate  the  over- 
lap bet-.y-een  x'  ^^-^  S.   The  process  is  then  repeated  for 
nev;  values  of  ^  until  this  overlap  is  a  maxiiuuin.   With 
this  final  ^  v.-e  set 

m  =  -2 


in  equation  (26)  and  solve  for  y 


CHF 


The  final  fit, 

therefore,  has  a  —■   v/eighting  factor,  introduced  to  remove 

r 
the  emphasis  that  the  r  =  3  Slater  orbital  inherently 

gives  to  the  tail  region.   In  practice  this  sometimes 
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serves  to  remove  slight  nodes  which  may  still  exist, 
(e.g.,  see  Figures  1,  2  and  3). 
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Figure  1.   Three  Different  Cu  4s  Orbitals. 

Curve  1  is  the  4s  orbital  from  the  full  SCF  calculation 
Curve  2  is  the  Coreless  Ilartree-Fock  orbital  with  the 
final  C  and  rn  =  0  (see  eqn ,  '  s  (25)  and  (26).   Curve  3 
is  the  final  CKF  orbital  (m  =  -2) . 
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B.   Generating  the  Coreloss  Hartree  Fock  Effective  Potential 

As  an  example  of  how  the  effective  potential  is  then 
obtained,  we  develop  here  the  Mn  's'  effe-  tive  potential. 
The  Hartree-Fock  Hamiltonian  for  the  4s  orbital  in  Mn  is 


"4s      2^ 


2, 

-  —  +  y  (2J  -K  )  +  j"^  +    y    (J- ,     -  ^K, ,    )    (28) 
r    ^     c   c      4s      i:  _    3d,m.    2  3d,m.       ' 


m^=-2 


which  we  replace  with 


H^HF  .   1,2  _  25  ^   FP 
4s      ^      r     4s 


m^=-2 


3d,mj^  -  2^3d,m^^   ^^^^ 


or 


4s      2      r     4s 


CFF  1 

in.=-2      'J.        '8, 


where 


FP  .  uEP  _  18 
4s     4s    r 


(31) 


HF       CHF 
Note  that  H.   and  H.    do  not  commute,  that  is  they  have 
4s       4s 

different  aigenfunctions  (for  the  "4s"  orbital  these  are 

HF        CHF 

0^   and  \.        respectively) ,  and  operate  on  slightly  dif- 
ferent spaces.   Equations  (20)  and  (30)  give 


.HF.  ,HF. 


HF.  HF. 


"4s'*4s>  =  ^^s\^^^ 


(32) 


and 


2? 


„CHF|,CHF^     CHFi.CHF^  ,^^. 

^A  \'t>  A  >     =     ^A  U'/i    ^  (33) 

4s'  4s       4s'4s 

EP       CHF 
At  this  point  V.   and  e.         are  not  known,  and  we  choose 
"^       4s       4s 


CHF  _   HF   _  ,-,.. 

e.        =   ^A         =    E.  (34) 

4s      4s      4s 


so  that  we  may  get  essentially  the  same  eigenvalue  spectrum 

^    „CHF     ^    „HF 
for  H     as  for  H 

EP 
V7e  nov7  wish  to  find  that  V.   which  will  generate 

CHF  .  .    . 

X-  "   and  e.  .   Thxs  s  potential  is  to  be  used  m  SCF  type 
4s        4s  ■' 

equations  where  its  effect  v/ill  be  felt  through  its  matrix 
elements.  Therefore  an  important  feature  of  the  effective 
potential  is  that  it  produce  the  proper  m.atrix  elements 

(rather  than  have  a  physically  appealing  form) .   Thus  the 

CHF 
best  way  to  determine  V     is  to  do  a  m.atrix  (in  this  case 

* 

vector) fit  over  some  arbitrary  basis  set.    Operating  on 

equation  (33)  from  the  left  v;ith  an  arbitrary  basis  set 
<  i I  we  get 


^.,„CHF,  CHF^    ^.|  CHF^  ,^^. 


which  combines  with  equation  (3  0)  to  give 

^   . |„EP|     CHF^  ,. I  ,    1„2.  7  ,^c^ 

<i|V43lx43    >    =    <iU43    -     (-7V    )    -    -  (36) 

4s  ^_„         3d,m  2    3d,m    ''^4s 

ill  «  ^  X/  X* 

* 

For  convenience,  the  arbitrary  basis  set  chosen  is  that 

^  ^    aHF        ^^   ^   CHF 
used  for  (J),   or  m  effect  y. 
^4s  '^4s 

2h 


C.   Functional  Form  of  the  Effective  Potential 


The  form  of  the  effective  potential  used  here  will 


be 


LMAX-1 
^        =  ^LMAX^""^    "^        I        I^V^(r)  l£m><)im|  (37) 


where 


AV^(r)    =  V^(r)    -   V^^^(r)  (38) 


and  v/here  V(r)     (=   V  (r)    or    AV    (r)  )    is   given   by 

2 

n.       -a.r 

V{r)    =    y   A^    r   "-    e      "■       .  (39) 

1 


The  parameters  of  equation  (39)  (A.,  n.  and  a-)  are  ad- 
justed to  give  the  best  solution  to  equation  (36)  (for  the 

EP 
r4n  s  potential)  for  V   by  performing  a  non-linear  least 

squares  fit.   To  see  how  this  form  is  arrived  at  we  note 
first  that  these  effective  potentials  have  an  essentially 
Gaussian  form  (equation  (39)).   This  will  allow  us  to  ar- 
rive at  simple  analytic  expressions  for  the  integrals  over 
Gaussian  basis  functions,  with  which  they  will  be  used 
(see  for  example  Reference  17or  Reference  29).   There  is 
also  an  I   dependence  which  is  not  cut  off  at  the  maximum 
i  value  of  the  core  orbitals.   This  satisfies  the  criterion 
described  earlier  for  the  £  dependence  of  an  effective 
potential.  We  do  note  however,  that  there  is  a  finite  cut 
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off  for  the  I   value.   This  can  best  be  explained  with  an 
example.   For  a  Mn  valence  orbital  we  expect  to  see  a 
very  different  potential  for  a  'd'  and  a  'p'  electron, 
where  for  Mn,  'd'  core  orbitals  do  not  exist.   For  an  'f 
orbital  we  expect  to  see  almost  the  same  potential  as  for 
'd'  since  neither  of  these  have  any  core  orbitals  of  the 
same  symmetry,  and  both  have  very  little  amplitude  in  the 
core  region  which  results  in  very  little  core  dependence 
of  the  coulomb  and  exchange  integrals  that  are  being  re- 
placed.  In  general,  it  should  suffice  to  set  LMAX  as  one 
greater  than  the  highest  i   value  of  the  core  orbitals. 
Thus  using  equation  (37)  for  Mn ,  orbitals  with  s,  p  and 
d  symmetry  on  the  effective  potential  center  will  all  feel 
different  potentials,  whereas  all  orbitals  v/hich  have  a 
symmetry  dependance  on  this  center  of  higher  Z   value  will 
feel  the  same  potential  as  d. 
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IV.   RESULTS 

Using  the  procedure  described  above,  Coreless  Hartree- 
Fock  Effective  Potentials  (CHFEP)  were  loveloped  for  atoms 
A   the  first  three  rows  of  the  periodic  table. 

A.   First  Row 

Tables  2  through  4  contain  the  CHFEP 's  for  Li 
through  Ne.   Each  of  these  contains  a  -n  /r  term  in  the 
d  part  of  the  potential  (that  felt  by  all  electrons) , 
which  just  represents  the  nuclear  attraction  term  of  an 
atom  with 


total     core     eff     v 


electrons.   In  addition  for  C,  U,    0  and  Ne  we  have  the 
two  terms: 

2  2 

T  ,   -80. Or      „  .  -30. Or 
-1. 6  e  -0. 4e 

r  '     r 

These  terms  v;ere  preset  at  these  values.   Together  they 
go  to  -2/r  (=  -n  /r)  as  r  goes  to  zero,  so  that  together 
with  the  -n  /r  term  above  there  is  an  overall  -z/r  depen- 
dence of  the  potential  as  r  goes  to  zero.   This  of  course 
introduces  a  constraint,  .;;...  :  „>  test  its  importance,  the 
flourine  CHFEP  was  generated  both  with  and  without  this 
constraint.   The  difference  in  the  two  effective  potential 
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was  found  to  be  negligible.   For  the  remaining  atons  it 
was  decided  to  include  this  constraint  for  aesthetic 
reasons  as  long  as  there  was  no  serious  difficulty  fitting 
a  potential  subject  to  this  constraint.   The  value  of  n 
as  well  as  the  number  of  terms,  for  the  different  poten- 
tials are,  as  described  above,  parameters  of  the  potential. 
Where  possible  these  were  kept  fairly  uniform.   Thus  for 
fcxam-ple  C,  N  and  0  are  alike  in  number  of  tarm.s  and  n 
values.   For  B  (v;hich  v/as  done  later)  it  v;as  found  neces- 
sary to  add  1  term  and  change  the  n  values  in  the  s-d  part 
of  the  potential  to  get  a  reasonable  fit.   For  Be  one  term 
for  the  p-d  potential  gave  essentially  as  good  a  fit  as 
tv;o  (this  is  primarily  because  the  p  and  d  electrons  see 
essentially  the  Scune  core  so  that  V  ^^^,  representing  the 
difference  between  the  p  and  d  potential  is  small) . 

To  see  how  well  these  potentials  perform,  the  ground 

state  and  various  excited  states  were  calculated  both  with 

21 
and  without  these  CHFEP ' s   .   The  results  of  these  tests 

are  presented  in  Tables  5  through  10 .   The  basis  sets  used 

to  test  these  were  the  IIS,  7P  basis  of  S.  Huzinaga  and 

y.  Sakai"'"^  contracted  to  (5,1,1,1,1,1,1/3,1,1,1,1).   In  all 

cases  the  errors  in  the  orbital  enr  .  lies  and  total  energies 

were  less  than  a  few  hundredths  and  a  few  thousandths  of 

an  a.u.  respectively  (or  a  fev;  tenths  and  a  few  hundredths 

of  an  electron  volt  respectively) .   In  most  cases  however, 

the  errors  v/ere  an  order  of  magnitude  better.   As  expected 

the  greatest  errors  occur  for  those  states  which  are 
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furthest  from  the  state  from  vrhich  the  potential  was 
created.   Thus  the  errors  in  the  energies  for  promoting 
an  electron  from  a  2s  orbital  to  a  2p  orbital  is  greater 
for  C  than  for  Be.   This  is  because  th  p  part  of  the 
potential  used  for  th~  carbon  configuration  (2s  2p  2p^2p  ) , 

V  :S   generated  from  a  p  orbital  in  the  ground  state 

2   11 
(2s  2p  2p  ) ,  whereas  for  Be,  in  order  to  get  the  p  part 

of  the  potential  the  p  orbital  of  the  2s  2p   configuration 
was  used,  so  that  it  v/as  generated  from  the  identical  state 
we  are  looking  at  (2s  2p  ) .   In  going  from  C  to  F  the  errors 
in  this  type  of  excited  state  are  comparable.   They  will 
however  be  most  serious  for  B,  C  and  ^7,  where  these  excited 
states  are  most  important  in  forming  molecules.   These  er- 
rors are  however  small  enough  that  no  serious  problems  are 
expected. 

To  see  hov;  well  the  tail  of  the  wave  function 

HF        CHF 
are  reproduced.  Figure  4  gives  ({>_   and  Xy        ^or   flourine. 

We  see  that  the  wavefunctions  too  are  well  described  in 
the  present  scheme. 

Finally  we  present  in  Table  11  a  list  of  4  term  s 
Gaussians  for  use  with  these  CIIFEP's.   V^ith  these  new  basis 
sets  the  number  of  functions  necessary  to  perform  the  SCF 
calculations  has  been  greatly  reduced,  resulting  in  major 
savings  in  computational  effort.   The  functions  given  in 
Table  11  were  obtained  by  non-linear-least-squarcs  fitting 
four  single  Gaussians  to  X-   and  keeping  the  resulting  ex- 
ponents.  For  F  and  0  these  basis  sets  were  then  tested, 
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and  the  results  are  given  in  Table  12.   The  basis  set 
used  for  these  tests  was  the  four  single  s  Gaussians  un- 
contracted  and  the  7P  Gaussians  of  S.  Huzinaga  and  Y. 
Sakai^^  contracted  to  (31111).  .The  results 
are  seen  to  be  quite  good,   demonstrating  that  the  pre- 
sent method  does  in  fact  accomplish  its  goal:  simplifica- 
tion of  SCF  calculations  without  introduction  of  serious 
errors. 
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TABLE  5.   Be  Tests^l^  Energies  in  Hartrees. 


State/Coupl: 

Lng 

Full  SCF 
-  0.30926 

Eff.  Pot 
-0.30931 

Error 

Be2s^ 

^2s 

.00006 

Singlet 

-14.57203 

-0.96772 

— 

Be2s^2pl 

^2s 

-  0.39369 

-0.39914 

.005 

Triplet 

"■2^ 

-  0.24138 

-0.24039 

.001 

E^ 

-14.51114 

-0.90606 

— 

AE^ 

0.06169 

0.06167 

.00002 

Be2p^2p^ 

^2p 

^T 

-  0.22060 

-0.22130 

.0007 

Triplet 

-14.33556 

-0.72542 

— 

AE^ 

0.23727 

0.24230 

.005 

Be2p-'-2p- 

'2p 

-  0.19074 

-0.19137 

.0006 

Singlet 

^T 

-14.30149 

-0.69131 

— 

AE^ 

0.27134 

0.27641 

.005 

Be2s-'-2p-'- 

^2s 

-  0.34220 

-0.35595 

.01 

Triplet 

^2p 

^T 
AE^ 

-  0.10971 

-0.10794 

.002 

-14.38617 

-0.78069 

— 

0.18666 

0.18704 

.0004 

Be"2s 

^2s 

-  0.08390 

-0.08427 

.0004 

Doublet 

^2p 

^T 

0.02876 

0.03064 

.003 

-14.53025 

-0.92340 

— 

AE^ 

0.04258 

0.04433 

.002 

Be"^2s^ 

^2s 

-  0.66611 

-0.67309 

.006 

Doublet 

^T 

-14.27719 

-0.67309 

AE^ 

0.26564 

0.29463 

.001 

Be2p^ 

^2p 

-  0.17637 

-0.17696 

.0006 

Singlet 

E^ 

-14.28492 

-0.67472 

— 

AE^ 

0.28791 

0.29300 

.005 
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TABLE  6.   B  Tests   ,  Energies  m  Hartrees, 


State/Coupling 


2   1 
B  2s  2p 

Doublet 


B  2s^2p^2pl 


Quartet 


B  2sl(a)2pl(6)2pl(3) 
Doublet 


B  2s^(B)2pl(a)2p^(3) 
Doublet 


B  2p^2p^2p^ 
Quartet 


B  2p^(a)2pl(n)2p^(B) 
Doublet 


B   2s^2p^2pl 


Full  SCF 

Eff.  Pot. 

Error 

^2s 

-  0.4947 

-0.49489 

.0002 

^2p 
^T 

-  0.30983 

-0.30921 

.0006 

-24.52871 

-2.55632 

— 

^2s 

-  0.64227 

-0.65079 

.007 

^2p 
^T 

-  0.35455 

-0.35637 

.002 

-24.45028 

-2.48427 

— 

AE^ 

0.07842 

0.07205 

.006 

^2s 

-  0.20659 

-0.19685 

.01 

^2p 

^T 

-  0.73404 

-0.74593 

.01 

-24.52587 

-2.54111 

— 

AE^ 

0.00284 

0.01521 

.01 

'2s 

-  1.08168 

-1.09611 

.001 

'2p(a) 

-  0.12583 

-0.12G03 

.0002 

^2p(B) 

-  0.88616 

-0.89233 

.003 

^T 

-24.66901 

-2.70195 

— 

AE^ 

-  0.14031 

-0.14563 

.005 

^2p 
^T 

-  0.32606 

-0.32906 

.003 

-24.12770 

-2.15507 

— 

AE^ 

0.40101 

0.40125 

.0002 

''2p(a) 
2p(B) 

-  0.80710 

-0.80683 

.0003 

-  O-l'-'f^ 

-0.16941 

.003 

^T 

-24.  .  ,  .; 

-2.51029 

— 

AE^ 

0.04473 

0.04603 

.002 

^2s 

-  0.22909 

-0.22S53 

.0006 

^2p 

^T 

-  0.02062 

-0.02030 

.0004 

-24.51665 

-2.54388 

-- 

AE^ 

0.01205 

0.01244 

.0004 
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TABLE  6   Continued 


State/Coupling 


Full  SCF   Eff.  Pot.   Error 


B   2s 


•2s 


Singlet 


T 


-  0.87382   -0.87545   .003 
-24.23722   -2.26568 
0.29148     0.29064    .0009 
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TABLE  7.   C  Tests   ,  Energies  in  Hartrees. 

State/Coupling  Full  SCF    Eff.  Pot.    Error 

2   11 
C  2s  2p  2p-^ 

■"  riplet 


2   11 
C  2s  2p  2p^ 


Singlet 


111" 
C  2s  2p  2p  2p" 

Quintet 


2s^{a)2p^(6)2p^(3)2p^(B) 
Triplet 


2s^(a)2p^(a)2p^(a)2p^(B) 
Triplet 


C  2s^(8)2p^(a)2p^(u)2p^{0) 
Singlet 


^2s 

-  0  '0560 

-0.70673 

.001 

^2p 
^T 

-  0.43332 
-37.68802 

-0.43269 
-5.35552 

.0009 

^2s 

-  0.71860 

-0.71986 

.001 

^2p 

-  0.38129 

-0.38062 

.0007 

-37.63073 

-5.29824 

— 

AE^ 

0.05729 

0.05723 

.0000 

^2s 

-  0.94171 

-0.95765 

.02 

^2p 
^T 

-  0.4788fi 

-0.48102 

.002 

-37.59860 

-5.27826 

— 

AE,j 

0.07941 

0.07726 

.002 

^2s 

-  0.19114 

-0.17366 

.01 

^2p 
^T 

-  1.63250 

-1.65272 

.02 

-38.75913 

-6.40821 

— 

AE^ 

1.07101 

1.05269 

.02 

^2s 

-  2.02937 

-2.04927 

.02 

''2p{a) 
^2p(0) 

-  1.75350 

-  0.11999 

-1.75835 
-0.12063 

.005 
.0006 

-38.98594 

-6.66397 

-- 

AE^ 

1.29792 

1.30346 

.01 

^28 

-  0.91148 

-0.93602 

.02 

*^2p(a) 
^2p(p) 
^T 

-  0.65796 

-  0.00472 

-0.63361 
-0.84191 

.02 
.04 

-37.97534 

-5.62634 

-- 

AE^ 

0.29733 

0.27082 

.03 

TABLE  7   Continuad 


State/Coupl 

ing 

Full  SCF 

Eff.  Pot. 

Error 

C  2s^2p'' 

^2s 

-  0.72534 

-0.72669 

.001 

Singlet 

^2p 

^T 

-  0.35659 

-0.35590 

.0007 

-37.60291 

-5.27043 

— 

AE^ 

0.08510 

0.08509 

.0001 

C~  2s^2p-'-2p 

lo  1 
2p 

^2s 

-  0.36466 

-0.36448 

.0002 

Quartet 

^T 

-37.70653 

-5.37368 

— 

AE^ 

0.01851 

0.01816 

.0003 

C^  2s22pl 

^2s 

-  1.15377 

-1.15711 

.003 

Doublet 

^2p 

^T 

-  0.90477 

-0.90311 

.002 

-37.29161 

-4.96001 

— 

AE^ 

0.39641 

0.39551 

.0000 
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TABLE  8.   N  Tests   ,  Energies  in  Hartrees. 


State/Coupling 


N  2s^2p^ 
Quartet 


12   11 
N  2s  2p  2p  2p 

Quartet 


N  2s^2p(a)2p(e)2p(B) 
Doublet 


2   2   1 
N  2s  2p''2p-^ 

Doublet 


2   2   11 
N   2s  2p  2p  2p 

Triplet 


Full  SCF 

Eff.  Pot. 

Error 

^2s 

-  0.941  2 

-  0.94659 

.005 

^2p 

-  0.56877 

-  0.56701 

.002 

-54.20416 

-  9.70829 

— 

^2s  . 

-  1.18936 

-  1.21378 

.02 

w 

-  0.42856 

-  0.42680 

.002 

^2p^ 

-  0.62783 

-  0.63311 

.005 

^T 

-53.80093 

-  9.31297 

— 

AE^ 

0.40322 

0.39533 

.008 

"2s 

-  0.90557 

-  0.91092 

.005 

^2d(8) 

-  1.20158 

-  1.19858 

.003 

^2p(a) 

-  0.34477 

-  0.34341 

.001 

^T 

-54.71478 

-10.21801 

— 

AE^ 

0.51062 

0. 50972 

.0009 

^2s 

-  0.96604 

-  0.97167 

.006 

%P^ 

-  0.55789 

-  0.55609 

.002 

^T 

-54.06590 

-  9.57019 

— 

AE^ 

0.13825 

0.13810 

.0002 

^2s 

-  0.56260 

-  0.56669 

.004 

^   2 

-  0  .  ^  -  ^  -^  -^ 

-  0.03031 

.00009 

2V 

%P^ 

-  0.14140 

-  0.14040 

.001 

=T 

-54.12460 

-  9.62790 

— 

«T 

0.07955 

0.08038 

.0008 

TABLE  8   Continued 


State/Coupling 

n"^  2s^2p^2p^    ^   2    "  O-^'^^O?    -0.94861    .007 

2p 
Triplet      e       -  1.07228    -1.07899    .007 
2p^ 
E       -52.67811    -8.16585 

AE„       1.52605     1.54244    .02 
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TABLE  9.   0  Tests   ,  Energies  in  Hartrees. 


State/Coupling 

2   2   11 
0  2s  2p  2p  2p 

Triplet 


"2s 
^2p^ 


Full  SCF 

-  1.24483 

-  0.56519 

-  0.70258 
-74.81001 


Eff.  Pot. 

-  1.24565 

-  0.56471 

-  0.70205 
-15.75898 


Error 

.0008 
.0005 

.0005 


12   2   1 
O  2s  2p  2p  2p 

Triplet 


2   2   11 
0  2s  2p''2p-^2p-^ 

Singlet 


2s 
^2p2 

^2pl 

^T 
AE„ 


2s 

^2p^ 

^T 
AE„ 


1.44663 
0.46232 


-74.18610 
0.62383 

-  1.25631 

-  0.57830 

-  0.62321 

-74.72811 
0.08189 


-  1.46498 

-  0.56336 


-  0.76812    -  0.77683 


-15.14453 
0.61440 

-  1.25727 

-  0.57777 

-  0.62266 

-15.67708 
0.08189 


.02 
.0006 

.009 


.009 

.0009 
.0005 

.001 
.000005 


2   2   2 
O  2s  2p  2p'' 

Singlet 


2s 

-2p 


AE, 


•  1.26232 

■  0.58518 

■74.63821 

0.12180 


-  1.26334 

-  0.58462 
-15.63717 

0.12181 


.001 
.005 

.000008 


2   2   2   1 
0   2s  2p  2p'^2p-^ 

Doublet 


^-2s 
2p 


■2p- 


*"T 
AE„ 


0.79822 
0.09676 


•74.78535 
0.02466 


-  0.79725 

-  0.09642 


-  0.20342    -  0.20307 


-15.73399 
0.02499 


.001 
.0003 

.0003 


.0003 


i.^ 


Table  9.   Continued 


State/Couplinq 

Full  SCF 

Eff.  Pot. 

Error 

-   2   2   2   2 
0  2s  2p  2p  2p 

^2s 

-  0.92771 

-  0.93744 

.01 

Doublet 

^2p 
^T 

-  0.11167 

-  0.11254 

.0009 

-74.10688 

-15.06066 

-- 

AE^ 

0.70313 

.69832 

.005 

■?-   2   2   2   2 
0   2s  2p'^2p^2p 

^2s 

-  0.44233 

-  0.44056 

.002 

Singlet 

'2p 

0.27485 

0.27519 

.003 

E^ 

-74.40161 

-15.34991 

— 

AE^ 

0.40839 

0.40907 

0"*"2s^2p-'-2p-'-2p-'- 

^25 

-  1.80072 

-  1.80449 

.004 

Quartet 

^2p 

-  1.32689 

-  1.32578 

.001 

1. 

^1 

-74.37125 

-15.32075 

— 

AE^ 

0.43875 

0.43823 

.0005 

hh 
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TABLE  10.   F  Tests   ,  Energies  m  Hartrees, 


State/Coupling 

Full  SCF 

Eff.  Pot. 

Error 

?   ?   2   1 
F   2s  sp  2p^2p-^ 

^2s 

-  1,57259 

-  1.57138 

.001 

Doublet 

%P^ 

%P^ 
^T 

-  0.70486 

-  0.70450 

.0004 

-  0.83442 

-  0.83409 

.0003 

-99.40858 

-24.00570 

— 

12   2   2 
F  2s-^2p^2p^2p^ 

^2s 

-  1.70537 

-  1.71821 

.01 

Doublet 

^2p 
^T 

-  0.70232 

-  0.70451 

.002 

-98.52932 

-23.13266 

— 

AE^ 

0.87926 

0.87304 

.006 

2   2   2   2 
F   2s  2p'^2p^2p^ 

"2s 

-  1.05757 

-  1.05420 

.003 

Singlet 

'2p 
-T 

-  0.16622 

-  0.16586 

.0004 

-99.45223 

-24.04S86 

— 

AE^ 

-  0.04365 

-  0.04317 

.0005 

f"*"  2s^2p^2p^2p^ 

"2s 

-  2.19507 

-  2.19730 

.002 

Triplet 

%P^ 

%P^ 
^T 

-  1.36156 

-  1.36107 

.0005 

-  1.52948 

-  1.52901 

.0005 

-98.83139 

-23.42901 

— 

AE^ 

0.57719 

0.57669 

.0005 
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TABLE  11.   Exponents  for  Four  Term  Gaussian  Basis  2s 

Coreless  Hartree  Fock  Orbitals.   Li  through  Ne. 
Li  Be  B  C 


2.89379 

1.92463 

3.48924 

5.43368 

0.61815 

C. 17264 

0.32044 

0.53387 

0.07385 

( .08894 

0.11488 

0.20524 

0.02817 

0.04000 

0.07004 

0.07958 

N  O  F  Ne 

7.81931  10.67015  13.94622  17.64863 

0.77843  2.06640  1.39517  1.77051 

0.29^76  0.41239  0.54084  0.08912 

0.11610  0.16019  0.20987  0.26649 
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TABLE  12.   Small  Basis  Set  Calculations  for  O  and  ^ ,    Energies 
in  Hartrees. 

+ 
State/Coupling         Large  Basis  Snail  Basis  Error 

F  2s^2p^2p^2p^ 
Doublet 


2   2   11 
0  2s  2p^2p-^2p-^ 

Triplet 


e_ 

-  1.57138 

-  1.55748 

.01 

2s 

e 

-  0.70450 

-  0.70720 

.003 

2p 

\^ 

-  0.83409 

-  0.83691 

.003 

^T 

-24.00570 

-23.9775 

.03 

^2s 

-  1.24565 

-  1.23333 

.01 

-  0.56471 

-  0.56762 

.003 

%  7 

-  0.70205 

-  0.70514 

.003 

2p 

^T 

-15.75898 

-15.73372 

.03 

Using  the  4s  type  functions  of  Table  II,  and  the  7p  type 
function  of  S.  Huzinaga  and  Y.  Sakai  (J.  Chen.  Phys.  50, 
1371  (1969))  contracted  (3,1,1,1,1). 
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B.   Second  Row 

The  effective  potentials  for  the  second  row  are  given 
in  Tables  13  through  15.   Analogous  to  the  CHFEP ' s  for 
the  first  row,  most  of  these  contain  the  term 


-10   -80. Or^ 
—  e 
r 


for  the  -n  /r  dependance.   Tables  16  -  21  contain  the  re- 
sults of  the  tests  for  this  row.   The  results  are  comparable 
to  those  of  the  first  row.   Since  d  type  polarization  func- 
tions are  important  for  the  second  row,  we  also  present 
some  tests  for  states  v/here  there  is  an  occupied  d  orbital. 
Figures  6a  and  6b  show  how  the  CHF  3s  and  3p  orbitals  for  Si  and 
AJl  compare  to  the  Hartree  Fock  3s  and  3p  orbitals.   The  s 

and  p  basis  sets  used  for  these  tests  was  the  12s,  9p  sets 

23 

of  Dunning   contracted  (2,1,1,1,1,1,1,1,1,1,1/ 

2,1,1,1,1,1,1,1) . 

Also  given  in  Table  22  are  the  basis  sets  for 
s  and  p  orbitals  for  use  with  effective  potentials.   Table 
23  gives  the  results  of  tests  v;ith  these  small  basis  sets 
for  S  and  CJ-. 
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Figure    7.       Sulfur   CHFEP ' s    for   S,    P   and   D. 

Curves 
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Curves  1,  2  and  3  are  Vl    ,    W  and  W,  respectively 
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21 
TABLE  16.   Mq  Tests   ,  Energies  m  Hartrees. 


S  tate/Coiipl  ing 

Full  5CF 

E^^.   T>ot. 

Error 

9 

Mg   3s" 

"3s 

-   0.25301 

-0.25301 

.000005 

Singlet 

^T 

-199.60837 

-0.78939 

—  — 

Mg   33-'-3p-'- 

"3s 

-   0.33991 

-0.34236 

.002 

Triplet 

^3p 
^T 

-   0.17939 

-0.17940 

.OOOni 

-199.54013 

-0.721f'9 

— 

AE^ 

0.06R24 

0.06740 

.0008 

Mg"^  3s  " 

"3s 

-   0.54145 

-0.54747 

.006 

Doublet 

^T 

-199.36545 

-0.54747 

— 

AE^ 

0.24292 

0.24192 

.001 

Mg   3s^3p-'-        c  -   0.10587  -0.106S5    .0007 

Doublet        e^  +   0.02276  +0.02359    .0008 

3p 

E  -199.58372  -0.76396 

AE^  0.02465  0.02543    .0003 

* 

Mg  3s-'"3d-'-         G^  -   0.^1993  -0.42617    .006 
^                   3s 

Triplet        £3^  -   0.05861  -0.05854    .0007 

E^  -199.24524  -0.60593 

AE^  0.36312  0. 183^^6    .2 


*This  state  v/as  done  v;ith  a  less  flexible  basis  set,  and 
comparisons  of  excitation  energies  are  not  strictly  valid, 
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TABLE  17.   AR  Tests^^.  Energies  in  Hartrees. 


Stato/Couplina 

Full  SCr 

E^.f.    Pot 

Error 

A?.  3s^3p-'- 

^3s 

-   0.39311 

-0.39254 

.0006 

Doublet 

"2p 

-   0.20982 

-0.21078 

.001 

^T 

-241.07124 

-1.89695 

— 

Al    3s-'-3p-'-3p-'- 

^3s 

-   0.51931 

-0.52255 

.003 

Quartet 

^2p 
^T 

-   0.26011 
-241.78531 

-0.26139 
-1.81594 

.001 

AE^ 

0.08593 

0.08101 

.005 

AZ'^   3s^ 
Singlet 

^35 

'^T 

-   0.C5222 
-241.66905 

-0.65592 
-1.69448 

.004 

AE^ 

0.20219 

0.20247 

.0003 

a;,"  3s^3p-'-3p-'- 

^3s 

-   0.18461 

-0.18350 

.001 

Triplet 

^3p 

-   0.00810 

-0.00719 

.0009 

^T 

-241.86649 

-1.89044 

-- 

AE^ 

0.00475 

0.00650 

.002 

2   1* 
At  3s  3d^ 

"33 

-   0.53057 

-0.53413 

.004 

Doublet 

"3d 

-   0.05804 

-0.05806 

.00001 

^•T 

-241.72705 

-1.75251 

— 

AE^ 

0.14419 

.14443 

.0002 

*See  footnote  of  Table  16. 
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21         .    . 
TABLE  18.   Si  Tests   ,  Eneraies  in  Hartrees. 


State/Couplinq 

Full  SCF 

Eff.  Pot. 

Error 

Si  3s^3p-^3p-^ 

^3s 

-   0.53931 

-0.53855 

.0008 

Triplet 

^3p 

-   0.29672 

-0.29791 

.001 

^T 

-288.84BG5 

-3.72568 

Si  3s-'-3p-'-3p-'-3p-'- 

^3s 

-   0.72274 

-0.72858 

.006 

Quintet 

^3p 

^T 

-   0.34823 
-288.75710 

-0.34948 
-3.64229 

.001 

AE^ 

0.09155 

0.08339 

.008 

Si"^  3s^3p^ 

^3s 

-   0.83851 

-0.84257 

.004 

Doublet 

^3o 

Ft 

-   0.58617 
-28  8.56722 

-0.59277 
-3,44441 

.007 

AE^ 

0.2814  3 

0.28127 

.0001 

Si~3s^3p^3p-'-3p^ 

^3s 

-   0.28322 

-0.23111 

.002 

Quartet 

'2p 

-   0.04916 

-0.04788 

.001 

^T 

-288.87772 

-3.75195 

— 

AE^ 

0.029075 

0.02628 

.003 

9   1   1  * 
Si  3s"3p-^3d-^ 

^3s 

-   0.70551 

-0.70972 

.004 

Triplet 

'3p 

-   0.45567 

-0.46223 

.007 

^3d 

-   0.06125 

-0.06120 

.00005 

^T 

-288.62835 

-3.50549 

— 

AE^ 

0.22030 

0.22018 

.0001 

*See  footnote  of  Table  16. 
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21 

TABLE  19.  ^P  Tests   ,  Energies  in  Hartrees. 


State/Coupling 

Full  SCF 

Eff.  Pot. 

Error 

„  -  2^  1-  1-  1 
P  3s  3p  3p  3p 

"3s 

-   0.69511 

-0.69522 

.0001 

Quartet 

'3p 

-   0.39088 

-0.39127 

.0004 

"  IT 

^T 

-340.64863 

-6,45646 

— 

P  3s  3p  3p-'3p-^ 

"3s 

-   0.87855 

-0.88937 

.01 

Quartet 

^T 

-   0.31316 

-0.31064 

.002 

-   0.44265 

-0.44635 

.004 

-340.35046 

-6.16225 

_— 

AE^? 

0.29817 

0.29421 

.004 

P"*"  3s^3p^3p^ 

"3s 

-   1.03358 

-1.03945 

.006 

Triplet 

^3p 

Em 

-   0.72032 

-0.72659 

.006 

-340.27978 

-6.08914 

-- 

AE^ 

0.36885 

0.36732 

.002 

„-  ,  2,  2-  1_  1 
P   3s  3p  3p  3p 

"3s 

-   0.41384 

-0.41353 

.0003 

%P^ 

%P^ 
^T 

-   0.02644 

-0.02324 

.003 

1^ 

-   0.09425 

-0.09292 

.001 

-340.62092 

-6.42278 

— 

AE^ 

0.02771 

0.03369 

.006 

2   111* 
P  3s  3p-^3p-^3d-^ 

"3s 

-   0.88540 

-0.89206 

.006 

^3p 
"3d 

-  0.57482 

-  0.06315 

-0.58153 
-0.06301 

.007 
.0001 

Km 

-340.33678 

-6.15178 

-- 

AE^ 

0.31186 

0.30460 

.007 

*See  footnote  of  Tabic  16, 
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TABLE  20.   S  Tests^^ 


,  Energies  in  Hartrees. 


State/Coupling 

2   2   11 
S  3s  3p  3p  3p 

Triplet 


F  3sl3p23p23pl 


'3  s 


3p' 
3p- 


'3s 


3p^ 


3p 


E 


T 

ae. 


Full  SCF 

-  0.87889 

-  0.39506 

-  0.47995 
-397.49870 

-  1.02846 

-  0.40457 


Eff.  Pot.    Error 


-397.07404 
0.42466 


■  0.87731 

■  0.39396 

•    0.48277 
•10.14770 

■  1.03734 

■  0.40380 


-   0.53564    -  0.54272 


9.72847 
0.41923 


.002 
.001 

,003 


.009 
.0008 

.007 


.005 


S"*"  3s^3p-'-3p-'-3p-'- 
Quartet 


'3s 

'3p 


AE. 


-  1.24043 

-  0.86231 
-397.16563 

0.33307 


1.24418  .0001 

0.87228  .01 
9.82035 

0.32735  .006 


2   2   2   1 
S   3s  3p  3p  3p 

Doublet 


'3s 


3p' 

3p 


AE. 


-  0.55877 

-  0.07730 

-  0.1^242 

-397.52367 
0.02497 


-  0.55544 

-  0.0740S 

-  0.14179 

-10.16538 
0.01769 


.003 
.003 

.0006 


.007 


2   1   1   1   1 
S  3s  3p  3p  3p  3d 

Quintet 


3s 
'3p 
'3d 


T 
AE. 


-  1.07332 

-  0.69951 

-  0.06519 
-397.23034 

0.26836 


1.08024  .007 

0.71205  .01 

0.06460  .0006 
9.88452 

0.26315  .005 


'See  footnote  of  Table  16. 
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21 
TABLE  21,  Ci   Tests   ,  Energies  in  Hartrees. 


State/Coupling 

^„  ^  2^  2_  2_  1 
Ca    3s  3p  3p  3p 

Doublet 


Full  SCF 


-3s 


3p' 

-   1 
3p- 


Eff.  Pot. 


Error 


-  1.07195  -  1.07169  .0003 

-  0.49011  -  0.48954  .0006 

-  0.56893  -  0.57240  .004 
-459.47461  -15.05678 


Ci 

12   2   2 
3s-^3p'^3p  3p 

^3s 

-   1.17465 

-  1.18327 

.009 

Doublet 

^3p 
^T 

-   0.49983 

-  0.50036 

.0005 

-458.90887 

-14.49352 

— 

AE^ 

0.56584 

0.56327 

.003 

CP' 

-  ^  2,  2,  2-  2 

3s  3p  3p  3p 

^3s 

-   0.71349 

-  0.71031 

.003 

Singlet 

'3p 

-   0.13271 

-  0.12919 

.004 

IT 

^T 

-459.56088 

-15.13495 

— 

AE^ 

0.08628 

0.07816 

.008 

CJ/ 

f  ,  2,  2-  1-  1 
3s  3p  3p  3p 

^3s 

-   1.47059 

-  1.47702 

.006 

Triplet 

%P^ 

%P^ 
^T 

-   0.88979 

-  0.89556 

.006 

-   0.99056 

-  1.00222 

.01 

-459.0/1087 

-14.62866 

— 

AE^ 

0.43384 

0.42812 

.006 

CJl 

3s^3p^3p-'-3p-'-3d-'- 

^3s 

-   1.29997 

-  1.30967 

.01 

Quartet 

%P^ 

-   0.72284 

-  0.73166 

.009 

-   0.82421 

-  0.83859 

.01 

-   0.06560 

-  0.06505 

.0005 

^'T 

-459.10603 

-14.69335 

— 

AE^ 

0.36857 

0.36344 

.005 

♦See  footnote  of  Table  16 . 
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TABLE  2 2 .   Exponents  for  Small  Gaussian  Basis  3s  and  3p 

Coreless  Hartree-Fock  Orbitals.   Na  through  Ar , 


Na 


Hq 


s 

P 

S 

P 

0.66980 
0.05367 
0.02222 

AJL 

1.48524 
0.07059 
0.02231 

1.10749 
0.08280 
0.01634 

Si 

3.13346 
0.14422 
0.03650 

%~f 

P 

S 

P 

1.53072 
0.13398 
0.01467 

P 

4.58868 
0.26631 
0.08118 
0.04543 

2.01139 
0.19160 
0.04268 

S 

6.17212 
0. 40375 
0.14448 
0.04925 

S 

P 

S 

-P 

2.48699 
0.26676 
0.08681 

Ci 

7.87999 
0.54427 
0.19410 
0.06181 

3.09265 
0.34139 
0.11756 

Ar 

9.44112 
0.81239 
0.35676 
0.12452 

S 

P 

S 

P 

3.75346 
0.41679 
0.14008 

11.67869 
0.91292 
0.33878 
0.12138 

4.29636 
0.52294 
0.18569 

13.45368 
1.18909 
0.48157 
0.16908 
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TABLE  23.   Small  Basis  Set  Calculations  for  S  and  CS,    ^   Enerai,  3 
in  Hartrees. 

State/Coupling  Large  Basis  Somali  Basis  Error 

S  3s^3p^3p-'-3p"'"    £3^  -  0.B7731  -  0.87602    .001 

Triplet        e   ^  -  0.39396  -  0.39334    .0006 

3p 

e   ,  -  0.48277  -  0.48172    .001 

3p^ 

E^  -10.14770  -10.14448    .003 

CH    3s^3p'''3p^3p-'-  c.^  -  1.07169  -  1.07171    .00002 

J  s 

e   „  -  0.48954  -  0.49014    .0006 

c  -  0.57240  -  0.57303    .0006 

3p 

E^  -15.05678  -15.05286    .004 

The  basis  sets  are  those  given  in  Table  22. 
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C,   Third  Row 

Tables  24-27  give  the  effective  potentials  for  K  through 

Zn.   Here  it  was  found  easier  to  fit  the  potentials  x>7ith- 

2 
out  the  (-n  /r)e      terra.   The  results  of  the  tests  given 

in  Tables  28  -  36  are  again  comparable  to  the  first  and 

second  rows.   The  basis  sets  used  for  these  tests  was 

24 
V7achters'    14s,  9p,  5d  contracted 

(6,1,1,1,1,1,1,1,1/5,1,1,1,1/4,1).   For  states  where  the 
4p  orbicals  were  occupied,  we  used  VJachters '  tv/o  extra  4p 
type  functions.   Again  we  present  basis  sets  for  s  func- 
tion in  this  rov;  in  Table  37.   These  have  been  tested  for 
Cu  and  Ti  and  the  results  of  these  tests  are  given  in 

Table  38 .   For  the  tests  of  these  functions  we  have  used 

24 
the  3-s  functions  given  in  Table  37  along  with  V7achters' 

5-d  functions  and  contracted  them  (1,1,1/5). 
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'D.OO 


Cu-Ws  ,  p,  D 


.80 


1.60       2.40 

R  n.u. 


3.20 


Figure  9.  "Copper  CHFEP ' s  for  5,  P  and  D. 

Curves  1,  2  and  3  are  W^ ,  Wp  and  V7^ 
respectively.   (VI^  =  ^"^g.^j  +  ^''^^' 
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TABLE  28.   Ca  Tests 


21 


Energies  in  Hartrees. 


State/Coupling 


Full  SCF 


Eff.  Pot.    Error 


Ca  4s 

Singlet 


■4s 


-   0.197 
-676.22511 


■0.196 
■0.61531 


.001 


Ca    4s-'"4p-'- 
Triplet 


4s 

^4p 

^T 

iE„ 


•  0.264 

•  0.150 
•676.18122 

0.04389 


-0.265 
-0.152 
-0.57542 
0.039G9 


.001 
.002 

.004 


Ca    4s    4p 

Sinqlet 


■4s 


T 
AE. 


-      0.237 
■      0.082 
■676.11694 
0.10817 


-0.244  .007 

■0.083  .001 
-0.509S6 

0.10545  .003 


Ca    4s-'-3d-'- 
Triplet 


'4  s 


■3d 


E 


T 
AE. 


■       0.182 
-      0.147 
-676.16243 
0.06263 


•0.185  .003 

-0.144  .003 
-0.55375 

0.06156  .001 


Ca    4s-'"3d-^ 
Singlet 


■4s 


3d 

T 

AE„ 


^T 


•      0.152 
■      0.147 
■676.14087 
0.08424 


■0.149 
•0.143 
■0.52653 
0.08878 


.003 
.004 

.005 
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TABLE  2  8   Continued 


State/Coupling 

Ca    4  s 

Doublet 


Full  SCP 

Eff.  Pot. 

Error 

^4s 

-   0.421 

-0.428 

.007 

^T 

-676.03/^^)3 

-0.42755 

— 

AE^ 

0.19018 

0.18776 

.002 

71 
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TABLE  29.   Sc  Tests   , 

Energies  in 

Hartrees. 

State/Coupling 

Full  SCF 

Eff.  Pot. 

Error 

Sc  4s^3d^ 

'4s 

-   0.208 

-0.208 

.000 

Doublet 

"3d 

-   0.34  4 

-0.342 

.002 

^T 

-759.70206 

-1.57242 

— 

Sc  4s-'-3d^ 

"4  s 

-   0.178 

-0.183 

.005 

Quartet 

"3d 

-   0.195 

-0.179 

.02 

^T 

-759.64580 

-1.51751 

— 

AE^ 

0.05626 

0.05491 

.001 

Sc  4s-'-4p-'-3d-'- 

"4s 

-   0.281 

-0.286 

.005 

Quartet 

"3d 

-   0.380 

-0.384 

.00'^ 

^4p 
^T 

-   0.158 
-759.66029 

-0.16  3 
-1.54019 

.005 

AE^ 

0.04177 

0.03223 

.01 

Sc-^^  4s4dl 

"4s 

-   0.448 

-0.460 

.01 

Triplet 

"3d 

-   0.600 

-0.606 

.006 

^T 

-759.50642 

-1.38186 

— 

AE^ 

0.19564 

0.19056 

.005 

Sc    4s 

"4s 

-   0.488 

-0.499 

.01 

Singlet 

^T 

-759.42857 

-1.28713 

— 

AE^ 

0.27349 

0.28529 

.01 

Sc-^1  3d2 

"4s 

-   0.442 

-0.431 

.01 

Triplet 

E^ 

X 

-759.46876 

-1.33549 

— 

AE^ 

0.23330 

0.23693 

.004 
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TABLE    2  9      Continued 


State/Coupl: 

Lng 

Full  SCF 

Eff.  Pot. 

Error 

Sc^2  3dl 

^3d 

-   0.904 

-0.923 

.02 

Doublet 

^T 

-758.06048 

-0.92314 

— 

AE^ 

0.64158 

0.64928 

.008 

^+2.1 
Sc    4s 

^4s 

-   0.762 

-0.793 

.03 

Doublet 

^T 

-758.94493 

-0.79320 

-- 

AE„ 

0.75713 

0.81529 

,06 

73 


21 
TABLE  30.   Ti  Tests   ,  Energies  m  Kartrees. 


State/Coupling 

Full  SCF 

Eff.  Pot. 

Error 

2   2 
Ti  4s  3d 

^4s 

-   0.218 

-0.219 

.001 

Triplet 

^3d 

-   0.438 

-0.434 

.004 

^T 

-848.36408 

-3.37018 

— 

9   2 
Ti  4s^3d^ 

^4s 

-   0.220 

-0.219 

.001 

Singlet 

^3d 

-   0.378 

-0.376 

.002 

^T 

-848.30326 

-3.30959 

-- 

AE^ 

0.06082 

0.06059 

.0002 

Ti  4s-'-3d^ 

^4s 

-   0.188 

-0.195 

.007 

Triplet 

^3d 

-   0.242 

-0.227 

.02 

^T 

-848.32167 

-3.33208 

— 

AE^ 

0.04  241 

0.03810 

.004 

„.  ,    1,    1^,2 
Ti  4s  4p  3d 

^4s 

-   0.298 

-0.309 

.01 

Triplet 

^3d 

-   0.473 

-0.482 

.01 

^T 

-   0.167 
-848.32477 

-0.165 
-3.33717 

.002 

AE^ 

0.03931 

0.03301 

.006 

Ti'^'-  4s^3d^ 

^4s 

-   0.473 

-0.490 

.02 

Doublet 

^3d 

-   0.705 

-0.707 

.002 

^T 

-848.16264 

-3.17730 

— 

AE^ 

0.20144 

0.19288 

.009 

Ti-'l  4s23d^ 

^4s 

-   0.512 

-0.522 

.01 

Doublet 

^3d 

-   0.900 

-0.9  33 

.03 

^T 

-848.01543 

-3.01288 

— 

AE^ 

0.34864 

0.35730 

.009 

7^ 


TABLE  3  0  Continued 


State/Coupling 

Full    SCF 

Eff.    Pot. 

Error 

Ti^l    3d^ 

^3d 

-       0.50^ 

-    .497 

.007 

Quartet 

^T 

-843.13518 

-3.13804 

-- 

AE^ 

0.22890 

0.23214 

.003 

m-+2       ot2 

Ti         3d 

^3d 

-      1.023 

-1.038 

.02 

Triplet 

^T 

-847.69154 

-2.68955 

— 

AE^ 

0.67254 

0.68063 

.009 

Ti^2    4313^1 

Us 

-      0.804 

-0.837 

.03 

Triplet 

^3d 

-      1.235 

-1.287 

.05 

^T 

-847.51270 

-2.50566 

— 

AE^ 

0.85138 

0.86452 

.01 
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21 
TABLE    31.      V  Tests       ,    Energies    in   Ilartrees. 


State/Coupling 

Full  SCF 

-Eff.  Pot. 

Error 

2   3 
V  4s^3d-^ 

"4s 

-   0.228 

-0.228 

.000 

Quartet 

"3d 

-   0.503 

-0.502 

.001 

^T 

-942.83485 

-6.14401 

— 

1   4 
V  4s-^3d 

"4s 

-   0.196 

-0.206 

.01 

Sextet 

"3d 

-   0.263 

-0.249 

.01 

^■T 

-942.80423 

-6.12117 

— 

AE^ 

0.0306 

0.0228 

.007 

V  3d^ 

"3d 

-   0.061 

-0.061 

.000 

Sextet 

^T 

-942.61667 

-5.92412 

— 

AE^ 

0.2182 

0.2207 

.003 

V  4s-'-4p-'-3d^ 

"4s 

-   0.314 

-0.327 

.01 

Sextet 

"3d 

-   C.546 

-0.546 

.000 

^4p 
^T 

~   0.171 

-0.169 

.002 

-942.79513 

-6.11378 

— 

AE^ 

0.0397 

0.0302 

.009 

V-*-^  4s23d2 

"4s 

-   0.534 

-0.543 

.009 

Triplet 

"3d 

-   1.020 

-1.050 

.03 

^T 

-942.45103 

-5.75215 

— 

AE^ 

0.3838 

0.3919 

.008 

V-^1  4sl3d3 

"4s 

-   0.491 

-0.517 

.03 

Quintet 

"3d 

-   0.781 

-0.784 

.003 

^T 

-942.62829 

-5.94  9  54 

— 

AE^ 

0.2066 

0.1945 

.01 
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TABLE  31   Continued 


State/Coupling 

Full  SCF 

Eff.  Pot. 

Error 

V-^l  3d* 

^3d 

-   0.537 

-0.533 

.004 

Quintet 

^T 

-942.60912 

-5.91671 

-- 

AE^ 

0.2257 

0.2281 

.002 

„+2  .  1^,2 
V    4s  3d 

^4s 

-   0.842 

-0.878 

.04 

Quartet 

^3d 

-   1.368 

-1.415 

.05 

^T 

-941.93104 

-5.23303 

— 

AE^ 

0.9038 

0.9110 

.007 

V^2  3^3 

^3d 

-   1.113 

-1.129 

.02 

Quartet 

^T 

-942.13388 

-5.43515 

— 

AE^ 

0.7010 

0.7096 

.009 
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21 
TABLE  32.   Cr  Tests   ,  Energies  in  Hartrees. 


State/Coupling 

Eull  SCF 

Eff.  Pot. 

Error 

Cr  4s-'-3d^ 

^4s 

0.204 

-  0.205 

.001 

Heptet 

^3d 

0.332 

-  0.312 

.02 

^T 

-1043.26492 

-10.01801 

— 

Cr  4s^3d'* 

^4s 

0.237 

-  0.228 

.009 

Quintet 

^3d 

0.537 

-  0.540 

.003 

^T 

-1043.24600 

-  9.99099 

-- 

AE^ 

0.01892 

0.02702 

.008 

Cr  3d^ 

"3d 

0.098 

-  0.078 

.02 

Quintet 

E^ 

-1042.91191 

-  9.652G0 

— 

AE^ 

0.35311 

0.36580 

.01 

Cr  4p-'-3d^ 

"3d 

0.398 

-  0.381 

.02 

Haptet 

"4p 

0.112 

-  0.113 

.001 

E^ 

-1043.17394 

-  9.9264 

— 

AE^ 

0.09098 

0.09200 

.001 

^  +1  .  2_,3 
Cr    4s  3d 

"4s 

0.554 

-  0.548 

.006 

Quartet 

"3d 

1.110 

-  1.144 

.03 

^T 

-1042.82924 

-  9.56355 

— 

AE^ 

0.43568 

0.45446 

.02 

Cr^l  Asha' 

"4s 

0.519 

-  0.52  5 

.006 

Sextet 

"3d 

0.825 

-  0.824 

.001 

^T 

-1043.03481 

-  9.79914 

-- 

AE^ 

0.23011 

0.21887 

.01 

78 


TABLE    32      Continued. 


State/Coupling 

Full  SCF 

Eff.  Pot. 

Error 

Cr    3d 

^3d 

0.618 

-0.603 

.02 

Sextet 

^T 

-1043.06212 

-9.81383 

— 

AF.^ 

0.20280 

0.20418 

.001 

„  +2  .  1^  ,3 
Cr   4s  3d 

^4s 

0.878 

-0.894 

.02 

Quintet 

^3d 

1.470 

-1.512 

.04 

^T 

-1042.29299 

-9.04440 

— 

AE^ 

0.97099 

0.97361 

.003 

cr*^  sa" 

^3d 

1.170 

-1.176 

.006 

Quintet" 

^T 

-1042.51832 

-9.27620 

— 

AE^ 

0.74657 

0.74223 

.004 
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TABLE  33.   Mn 

Tests   , 

Energies  in 

Hartrecs. 

State/Coupling 

Full  SCF 

V.^f.     TDot. 

Error 

Mn  4s^3d^ 

^4s 

0.245 

-  0.244 

.001 

Sextet 

^3d 

0.628 

-  0.627 

.001 

^T 

-1149.78740 

-15.21117 

— 

Mn  4s-'-3d^ 

^4s 

0.207 

-  0.215 

.008 

Sextet 

^3d 

0.357 

-  0.347 

.01 

^T 

-1149.62855 

-15.06271 

— 

AE^ 

0.15884 

0.1^846 

.01 

Mn  3d'^ 

"3d 

0.087 

-  0.085 

.002 

Quartet 

^T 

-1149.30734 

-14.73^80 

— 

AE^ 

0.48006 

0.47137 

.01 

Mn  4s  4p  3d 

"4s 

0.342 

-  0.361 

.02 

Octet 

^4p 
"3d 

0.179 

-  n.l76 

.003 

0.666 

-  0.669 

.003 

^T 

-1149.74638 

-15.18506 

— 

AE^ 

0.04102 

0.02611 

.02 

m^-^    4s-'-3d^ 

"4s 

0.539 

-  0.566 

.03 

Heptet 

"3d 

0.925 

-  0.927 

.002 

^T 

-1149.57207 

-15.01377 

— 

AE^ 

0.21533 

0.19740 

.02 

Mn-^1  3d^ 

"3d 

0.652 

-  0.652 

.001 

Quintet 

E-r 

-1149.42283 

-14.8^861 

— 

AE^ 

0.36457 

0.36256 

.002 

+  1    7   4 
Mn  -^  4s  3d^ 

"4s 

0.574 

-  0.582 

.008 

Quintet 

"3d 

1.216 

-  1.241 

.03 

^T 

-1149.29416 

-14.71072 

— 

AE^ 

0.49324 

0.50045 

.007 
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TABLE  3  3   Continued 

State/Coupling  ^ull  SCF  Eff.  Pot.    Error 

Mn"^-"-  3d^         3^  -    1.284  -  1.299      .02 

Sextet       j.^  -1149.03503  -14.45092 

E„                   0.75237  0.76025    .008 
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TABLE  34.   Co 

Test: 

21 

5   ,  Energies 

in  Hartrees. 

State/Couplinq 

Full  SCF 

Eff.  Pot, 

Erro: 

Co  4s^3d'^ 
Quartet 

^4s 
"3d 

0.263 
0.616 

-  0.263 

-  0.614 

.000 
.002 

^T 

-1381.29988 

-29.50895 

— 

Co  4s-'-3d^ 

"4s 

0.214 

-  0.218 

.004 

Quartet 

"3d 

0.348 

-  0.?42 

.006 

^T 

-1381.20349 

-29.41772 



AE, 

0.09649 

0.09123 

.005 

9 

Co  3d^ 

"3d 

0.068 

-  0.069 

.001 

Doublet 

^T 

-1380.88692 

-29.10359 

— 

AE^ 

0.41306 

0.40536 

.008 

Co  4p-'-3d^ 

^4p 

0.113 

-  0.114 

.001 

Quartet 

"3d 

0.426 

-  0.427 

.001 

^T 

-13B1.103C3 

-29.31405 

— 

AE^ 

0.19695 

0.19490 

.002 

Co-^1  ^shd' 

"4s 

0.568 

-  0.587 

.02 

Quintet 

"3d 

0.931 

-  0.932 

.001 

^T 

-1381.05908 

-29.27848 

.01 

Co-^1  3d^ 

"3d 

0.661 

-  0.664 

.003 

Triplet 

^T 

-1380.98999 

-29.20002 

— 

AE^ 

0.30999 

0.30893 

.  01 

^  +1  ,  2-, 6 
Co    4s  3d 

"4s 

0.613 

-  0.621 

.008 

Quintet 

"3d 

1.191 

-  1.206 

.02 

^T 

-1380.90307 

-29.10926 

— 

AE^ 

0.39691 

0.39969 

.000 
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TABLE  34   Continued 


State/Coupling  Full  SCF  Eff.  Pot.    Error 

Co"^^  Sd"^        e^^  -    1.316  -  1.330      .01 

Quartet      E^  -1380.49332  -23.69366 

AE^        0.80666      0.81529    .008 

Co"^-^  3d^        E^-  -    2.006  -  2.048      .04 

id 

Quintet      E^  -1379.34173  -27.51773 

AE^        1.95825      1.99122    .03 
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21 
TABLE    35.       Cu   Tests       ,    Fnerqies    m  llartrees. 


S tate/Coupl ing 

Full  SCF 

Eft:.  Pot. 

Error 

CU  Ashd^' 

^4s 

0.22918 

-  0.22957 

.0004 

Doublet 

^3d 

0.46044 

-  0.46145 

.0001 

^T 

-1638.80673 

-50.69719 

— 

CU  4pl3d^° 
Doublet 

sd  2 

0.11653 
0.54425 

-  0.11665 

-  0.54852 

.0001 
.004 

z 

-  T 

-1638.69498 

-50.58495 

-- 

AE^ 

0.11174 

0.11223 

.001 

2   9 
Cu  4  s  Sd"^ 

^4c 

0.28306 

-  0.23152 

.002 

Doublet 

%d^ 

%d^ 
^T 

0.83326 

0.75138 

-1638.80610 

-  0.84104 

-  0.75837 
-50.69201 

.009 
.007 

AE^ 

0.00027 

0.0052 

.0054 

CU-^^  4sl3d^ 
Triplet 

^4s 
%d^ 

%d^ 
^T 

0.59826 
1.09102 

1.17499 

-1638.53760 

-  0.60691 

-  1.10079 

-  1.18724 
-50.42724 

.009 
.01 

.01 

AE^ 

0.26912 

0.26995 

.0008 

cu+i  3a^° 

^3d 

0.78921 

-  0.79435 

.005 

Singlet 

^T 

-1638.57854 

-50.46837 

-- 

AF,^ 

0.22819 

0.22882 

.0006 

+  2    9 
Cu    3d-^ 

Doublet 

%dl 

%d^ 
^T 

1.57043 

1.45452 

-1637.94136 

-  1.5887 

-  1.4717 
-49.82229 

.02 
.02 

AE^ 

0.86536 

0.87490 

.01 
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21 
TABLE  36.   Zn  Tests   ,  Energies  in  Hartrees. 


State/Coupling 

^4s 

Full  SCF 
0.286 

Eff.  Pot. 
-  0.286 

Error 

zn  Asha'-'' 

.000 

Singlet 

^3d 

0.755 

-  0.754 

.001 

^T 

-1777.67477 

-63.89602 

— 

Zn  4s-'-4p-'-3d-'-° 

^4s 

0.386 

-  0.390 

.004 

Triplet 

^3d 

0.026 

-  0.826 

.000 

"4p 

0.185 

-  0.184 

.001 

^T 

-1777.57980 

-63.30177 

-- 

AE^ 

0.09497 

0.09425 

.0007 

„    +1  ,    1,,10 
Zn    4s  3d 

^4s 

0.604 

-  0.614 

.01 

Doublet 

^3d 

1.102 

-  1.105 

.003 

^'V 

-1777.39988 

-63.62244 

— 

AE^ 

0.27489 

0.27358 

.001 

^  +1  .  2_,9 
Zn    4s  3d 

^4s 

0.669 

-  0.676 

.007 

Doublet 

^3d 

1.416 

-  1.430 

.01 

^T 

-1777.12848 

-63.34657 

— 

AE^ 

0.54629 

0.54  94  5 

.003 

zn*'    3d" 

^3d 

1.520 

-  1.533 

.01 

Singlet 

^T 

-1776.79711 

-63.00994 

-- 

AE^ 

0.87766 

0.88608 

.008 
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TABLE  3  6   Continued. 


State/Coupling  Full  SCF  Eff.  Pot.  Error 

Zn"*"^  4s-^3d^      e^g  -    1.048  -  1.071  .02 

Triplet      £3^  -    1.843  -  1.865  .02 

E  -1776.47022  -62.68567 

AE^         1.20455       1.21035  .006 
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TABLE  37.   Exponents  for  Three  Term  Gaussian  Basis  4s 

Goreless  Hartree-Fock  Orbitals.   K  through  Zn. 


K  Ga  Sc         Ti 


0.270275 

0.44883 

0.63200 

0.6780 

0.040528 

0.06209 

0.03132 

0.0651 

0.018283 

0.02504 

0.02834 

0.0209 

V  Cr  Mn  Fe 

0.75580  0.76559  0.87877  0.97543 

0.07336  0.07017  0.08105  0.08824 

0.02590  0.02638  0.01339  0.01504 

CO  Ni  Cu  Zn 


1.08436 

1.1760 

1.0720 

1.25022 

0.09910 

0.1110 

0.0995 

0.13642 

0.03285 

0.0393 

0.0357 

0.04718 
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TABLE  38.   Small  Basis  Set  Calculations  for  Ti  and  Cu  ^  Energies 


in 

Hartrees. 

State/Coupl 

ing 

Larqe  Basis 

Small  Basis"*" 

Error 

2   2 

Ti  4s  3d 

Triplet 

^4s 

^3d 

-  0.218 

-  0.437 

-  0.219 

-  0.434 

.001 
.003 

^T 

-  3.36629 

-    3.36976 

.003 

CU  4sl3dl° 
Doublet 

^4s 
^3d 

-  0.22928 

-  0.45950 

-  0.22969 

-  0.46391 

.0004 
.005 

^T 

-50.69719 

-50.69714 

.00005 

"""using  the  3s  type  functions  of  Table  37  vzith  the  5d  type 
functions  of  V7achters  '   contracted  (1,1,1/5). 
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D.   Cu  One  Electron  CHFEP 

The  effective  potentials  we  have  developed  until 
this  point  were  based  on  the  usual  core-valence  separation 
of  orbitals.   We  know  however  that  in  some  cases,  orbitals 
that  are  usually  classified  as  valence  orbitals,  can  be 
considered  core  orbitals.   Mary-Beth  Baillie  et  al.   '   ' 
have  found  for  example  that  for  Cu  in  CU2,  CuH  and  Cu-^H^, 
the  d  orbitals  of  Cu  do  not  undergo  any  significant  change 
upon  molecular  formation.   It  would  then  be  an  additional 
savings  in  computational  effort  if  the  d  orbitals  v;ere  in- 
corporated into  the  effective  potential.   This  one  elec- 
tron potential  along  v/ith  some  tests  are  presonted  in 
Tables  3  9  and  4  0  respectively.   Inclur^ed  in  Table  39  is 
the  Ni  one  electron  effective  potential  developed  by  Carl  F. 
Melius.   Marv-Reth  Baillie,  et  al.^^'^^'^^  ^^^^  ^^^^^^  ^^^^ 
for  CU2,  CuH  and  CU2H2  this  one  electron  effective  potential 
reproduces  the  results  of  more  elaborate  calculations  quite 
well. 
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TABLE 

39. 

Ni  and  Cu  One 

Electron  1 

Ni 

n 

i 

A 

^S 

-2 
0 

1.17326 
0.12656 

2.60210 
-  0.06672 

Vd 

-2 

0.01299 

0.15429 

-1 

0.0 

-  1.0 

Vp-d 

-2 
0 

1.08491 
0.06158 

Cu 

1.79801 
-  0.01714 

n 

1 

A 

^s-d 

-2 

1.15858 

2.96371 

0 

0.25712 

-  0.19867 

Vd 

-2 

0.03588 

0.14945 

-1 

0.0 

-  1.0 

^-d 

-2 
0 

3.59990 
0.12943 

12.44478 
-  0.05452 

22 
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TABLE  40.   Cu  One  Electron  Effective  Potential  Tests,  Energies 
m  Hartrees. 

^4s        ^-^^"^    ^4p        ^^^Q^ 

Full  Basis 

without  Potential   -0.22940     —     -0.12087 

Full  Basis 

with  11  Electron 

Potential  -0.22357    .0003    -0.12049    .0004 

Full  Basis  v;ith 

1  Electron 

Potential  -0.22S66    .0007    -0.12053    .0003 

Three  Term  Basis 
v;ith  11  Electron 
Potential  -0.22969    .0003 

Three  Term  Basis 

v;ith  1  Electron 

Potential  -0.22850    .0009 


Using  the  basis  set  given  in  Table  37. 
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V.   CONCLUSION 

The  effective  potentials  that  have  been  presented  and 
tested  here,  appear  to  provide  a  promising  approach  to  SCF 
calculations.   The  essential  features  of  complete  SCF  cal- 
culations are  reproduced  quite  accurately.   The  agreement 
between  the  effective  potential  SCF  calculations  and  com- 
plete SCF  calculations  for  various  excited  states  show 
that  the  underlying  approximation  of  the  frozen  core  does 
not  introduce  any  serious  error. 

In  addition,  by  reducing  the  basis  sets  to  a  size 
needed  by  just  the  valence  pseudoorbitals,  v/e  are  able  to 
greatly  reduce  the  time  and  cost  of  SCF  calculations.   It 
appears  to  be  worthv/hile  to  continue  the  development  of 
CHFEP ' s  for  the  second  and  third  row  transition  elements. 
No  great  difficulty  is  forseen  in  doing  so. 
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•^x  "^y  '^  z 
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""  2p-^ 

the  singly  occupied  2p   and  2p   orbitals.   Similarly, 

for  the  transition  elements  many  of  the  d  orbitals 
with  different  spatial  s\'mm.etry  (dxy,  dxz ,  dyz,  d  -   -, 

X  -y 
d  2)  have  different  orbital  energies.   For  those  cases 
z 

v/e  present  only  one  each  of  a  singly  occupied,  and  a 
doubly  occupied  orbitals,  or  just  one  sample  d  orbital 
energy. 

94 


22.   The  effective  potentials  presented  in  these  tables 
(w   )  are  for  the  effective  Ilamiltonians 

H^P  =  T  +  T^^^^  +  Vy   =  T  -  -^  +  V^^  +  V^ 

A         r  X 

V  V 
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^v 
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LMAX-1 
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1     1 


VHr)    =    I  A^   r   ^   e 
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